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FOURIER-BESSEL HEAT KERNEL ESTIMATES 


JACEIC MALECKI, GRZEGORZ SERAFIN AND TOMASZ ZORAWIK 


Abstract. We provide sharp two-sided estimates of the Fourier-Bessel heat kernel and 
we give sharp two-sided estimates of the transition probability density for the Bessel 
process in (0,1) killed at 1 and killed or reflected at 0. 


1. Introduction 


We consider the Fourier-Bessel heat kernel, which is represented in terms of the Bessel 
functions of the first kind J v {z) and its successive n-th positive zeros \ n>v in the following 
way 


G u t 0, V ) = 2 (xy) u ex P ( _A n,^) 

71=1 


Jv Jij{^n,vy') 

\Ju + l{Ku )\ 2 


x,t/G(0,l), t > 0, (1.1) 


where v > — 1. The main results of the paper are the following sharp two-sided estimates 
of Gt(x,y ) given in 


Theorem 1. For every v > 
(1 + t) u+2 


Gt(x,y) 


1 we have 

{l-x)(l-y) 


1 A 


(t + xy) u+1 / 2 
whenever x, y G (0,1) and t > 0. 


t 


7t exp 


f \ x — y I 2 



( 1 . 2 ) 


V 

Here ~ means that the ratio of the functions on the right and left-hand side is bounded 
from below and above by positive constants depending only on v. Since the sum in (II.Ill 
is oscillating, this explicit representation can be only used to examine the behaviour of 
the kernel for large times. Indeed, it well-known that G v t (x, y) behaves like the first 
term of the series, whenever t > T 0 > 0 (for every T 0 , if we consider the upper-bounds 
and for some T 0 , if we deal with the lower-bounds). However, the description of the 
behaviour of G”(x,y) for small times is very difficult to obtain from the above-given 
series representation, since the sum is highly oscillating and the cancellations between 
the terms matter in that case. This is a reason why we do not use (11.111 in examine the 
small-time behaviour, instead we explore the relation between the Fourier-Bessel heat 
kernel and the transition probability density of the Bessel process with index v reflected 
at 0 and killed at 1. This approach enables us to use probabilistic tools like, for example, 
the Hunt formula or the Strong Markov property, but still the purely analytic studies of 
the properties of the modified Bessel functions are crucial for the proofs. 

The Fourier-Bessel expansions naturally associated with the Fourier-Bessel heat kernel 
G v t ( x , y ) has been studied for a long time in many different contexts, such as the study 
of the fundamental operators associated with the Fourier-Bessel expansions (see 0. H. 
0 , 0 , 0) or the related Hardy spaces (0) just to list a few from the latest works 
(see m for more references). Moreover, the Fourier-Bessel expansions are successfully 
applicable in variety of areas outside Mathematics. The estimates of G u t ( x , y ) has been 
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recently studied in m and pj2] , where the provided two-sided estimates of G" (x, y) were 
quantitatively sharp, i.e. the different constants appear in the exponential terms of the 
lower and upper bounds, ft makes the estimates not sharp, whenever \x — y | 2 >> t. In 
the estimates given in Theorem Q] the exponential behaviour of the kernel is described 
explicitly, i.e. the exponential terms in the lower and upper bounds are exactly the same. 
Such accurate results seem to be quite rare. Notice that even in the classical setting of 
Laplacian in R n , the known estimates of related Dirichlet heat kernel for smooth domains 
(see nan are also only quantitatively sharp (see also [13] and the references therein for 
corresponding results on manifolds). However, in the recent papers [I] and [2] the sharp 
two-sided estimates for the Dirichlet heat kernel of the half-line (a, oo) associated with 
the Bessel differential operator has been obtained. 

As we have previously mentioned, the result can be equivalently stated in the prob¬ 
abilistic context. More precisely, if we denote by p^\t,x,y) the transition probability 
density (with respect to the speed measure m^ u \dy) = y 2u+1 dy) of the Bessel process 
with index v > — 1 killed at 1 and reflected at 0, then we have p^\2t,x,y) = G^(x,y) 
and consequently 


Corollary 1. For given u > — 1 we have 


Pi\t,x,y) « 


(1 + t) v+2 
(t + xy) u+1 / 2 


q . t 1 - apt 1 -v) \ 


vt exp 


f \x — y\ 2 
(- 2t 


whenever x, y G (0,1) and t > 0. 



Futhermore, instead of studying the Bessel process reflected at 0, we can impose killing 
condition at both ends of the interval (0,1). Then we can expend the range of the 
index of the process to the whole real line (p. G R) and denote by ^(£, x, y) the 
transition probability density (with respect to the speed measure mS^\dy) = y 2fl+1 dy) of 
the corresponding process, i.e. the Bessel process killed when it leaves (0,1). Note that for 
p > 0 the process does not hit 0. Consequently the condition at zero (killing or reflecting) 
is relevant for the considered problem in that case, which means that p[oi)(^, x i u) and 

Pi'\t,x,y ) are identical for p = v > 0. Moreover, for p < 0 we can use the absolute 
continuity property of the Bessel process with different indices to get that for every p > 0 
we have 

p\p!i)( t > x ’y) = M 2 *PQ\){t,x,y), x,ye( 0,1), t> o. 

Collecting all together we obtain 


Corollary 2. For given p G R we have 



(1 + t) M+2 
(■t + xy)M +1 / 2 


\ A (l -s)(i -y) 

t 


1 A (xy) 2/i ( \x — y\ 2 


A j v t ^ 

(1.4) 


whenever x, y G (0,1) and t > 0. 


The paper is organized as follows. In Preliminaries we collect some basic properties 
of the modified Bessel function of the first kind I v (z) together with the estimates of the 
ratio of the form I u (y)/I u (x). Then we introduce the basic notation together with some 
properties of the Bessel processes, which are used in the sequel. In Section [3] we provide 
the proof of Theorem [T] which is divided into two parts. The first one relates to the lower 
bounds and the estimates in that case are given in Proposition [D The upper bounds are 
proved in Proposition [2j 
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2. Preliminaries 


2.1. Modified Bessel functions of the first kind. The modified Bessel function of 
the first kind is defined by 


(zy +2k i 

k\r(k + u+iy 


z > 0, v > —1. 


The above-given definition immediately implies that for every v > — 1 we have 

)■ *^ 0+ - P' 1 ) 

Moreover, the behaviour at infinity is described by 

Iv{z) -4=, 2 ->• OO. (2.2) 

\/ Zttz 

Here ~ means that the corresponding limit of the ratio of the both functions is 1. Finally, 
we will need the upper bounds of the ratio I v (y)/I v {x), which can be found in [9] 

9>z>0, ,>-1/2. (2.3) 

However, the above-given result is not true for v G (—1, —1/2). Since we will need such 
kind of estimates only for large x, y, we introduce the following result valid for every 
v > — 1 and large x and y. Note that the exponential term in the upper-bounds is the 
same as in (12.3(1 . but the factor ( y/x) u is here replaced by (y/x) v+1 . 


Lemma 1. For every y > x > 1 and v > — 1 we have 


< p y-x 

I v {x) ~ 



(2.4) 


Proof. As we have mentioned, it is enough to consider v G (—1, —1/2) and the proof in 
this case will be a slight modification of that given in |9]J. For every z > 0 we have (see 

cm) 

Iu+l(z) Z 

Iu{z) z + v + 1/2 
Thus, for every z > 1, we can write 

h±iW < —£— <i + l. 

I„(z) z — 1/2 z 

Consequently, the recurrent relation for the modified Bessel function implies 

zl' u {z) = vl v (z) + zl u+ i(z) < (u + 1 )I„{z) + zl v (z). 


Dividing both sided of the above-given inequality by zl u (z) and integrating the obtained 
relation lead to 



dz < {v + 1) 


f y dz 

/ T+y 


whenever y > x > 1. This gives 


IM < 

Iv(x) ~ 



e y ~ x , 


which ends the proof. 


□ 
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2.2. Bessel processes. We write Px and Ex'* for the probability law and the ex¬ 
pected value of a Bessel process with an index v e R on the canonical path space 
with starting point x > 0. The filtration of the coordinate process R(t) is denoted by 
T t = cr{i?(s) : s < t}. The transition density function (with respect to the speed mea¬ 
sure mS v \dy ) = y 2u+1 dy) of the process (with reflecting condition impose on 0, when 
v G (—1,0)) is given in term of the modified Bessel function in the following way 

{xy)~ 


P {u) (t,x,y) = 


■ exp 


x 2 + y 2 

2 T~ 


x, y > 0, t > 0. 


and 


p {v) (t,o,y) = 


exp 


xy_ 
t 

■ftI’ y - t>0 - 


(2.5) 


( 2 . 6 ) 


2 u t u+1 T(v + 1) 

Taking into account the asymptotic behavior of I v (z) at zero (12.11) and at inhnity (12.lj) 
we obtain 


p M (t,x,y) 


V 


: exp 


(x - yf 
2 1 


(2.7) 


(xy + t) u+l / 2 y/t 

whenever x, y > 0 and t > 0. 

Let us denote the first hitting time of a given point a > 0 by 

T a = inf(7 > 0 : R{t) = a} 

and we write qxl{ds) for the associated probability distribution with respect to Px *. 
Moreover, we introduce the first exit time from the interval (a, b ) 

T( a ,b) = inf {t > 0 : R(t) (a, &)}, 0 < a < b. (2.8) 

The laws of Bessel processes with different indices are absolutely continuous and the 
corresponding Radon-Nikodym derivative is described by 


dP^ 


dP 


iy) 


Tt 


m 


x 


(i—i/ 


exp 




ds 


R 2 (s ) 


(2.9) 


where x > 0, /x, v e R and the above given formula holds P^-a.s. on {To > t}. 

The transition probability density function p^\t, x, y) of the process killed at 1 can be 
expressed by the Hunt formula in the following way 

Pi ] (t, x, y) = p (u) (■ t , x, y) - r[ u) ( t,x,y ), (2.10) 

where 

rP(t,x,y) = EP[t<T 1 -,p^(t-T 1 ,R(T 1 ),y)]= f p^(t-s,l,y)&ds). 

Jo 

The last equality follows from the continuity of the path, i.e. the fact that R{T\) = 1 
Pi^-a.s.. 

Finally, we denote by P^ b ^{t,x,y) the transition probability density function of the 
process killed, when it leaves the interval (a, b), for given 0 < a < b, i.e. 

^ < mM(ck,) t) 6 xe ^ b ^ t>0 - 

We will denote the index of the process by /i, (/i 6 R), when we deal with the process 
killed at 0 (which is a case when a = 0 above) to distinguish this situation from the 
case when 0 is reflecting. Since the considered Bessel process with index —1/2 can be 
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represented as the norm of one-dimensional Brownian motion, p^ 0 ^ 2 \t,x,y) coincides 
with the corresponding object for Brownian motion on the real line, i.e. we have 


p{ 0 ^ 2) {i, x ,y) = 


\/ 2 iit _ 


E 


exp 


(x — y + 2 n,y 
2 1 


— exp 


(x — y — 2 n) 

2T 


2 \ 1 


(2.n) 

Moreover, the sharp two-sided estimates of p^^ 2 \t,x,y) are of the form (see Pyc, 
Serahn, Zak) 

{x-yf' 


p[ 0 ^ 2 \t,x,y) ~ ( 1 A 


xy 


(1 — t)(1 — y)\ 1 

1 A ----— -7= exp 

t J yft 


2 1 


( 2 . 12 ) 


whenever x, y G (0,1) and t < 1. Using the scaling property and the shift-invariance of 
one-dimensional Brownian motion, we arrive at 

n ( " 1/2) U x v) = -^—v { ~ 1/2) ( _-_ X ~ a y ~ a \ 

P(a,b) \hx,y) b _V ( o,i) _ a )2’ b _ a ’ b _ a ) ’ 


for given 0 < a < b. It leads to 


HJ) /2) m^) ~ ( 1 A 


(x - a)(y - a ) 


1 A 


(b-x)(b-y)\ J_ 

Vt 


exp 


{x-yf 

2 1 


(2.13) 


3. Proof of Theorem Q] 


The proof of Theorem Q] is divided into two parts, the first one relates to the lower 
bounds (Proposition [1]) and the other is devoted to show the upper bounds (Proposition 
ED. Moreover, due to the symmetry, we will generally assume in the proofs in this section 
that y > x. 

Proposition 1. For every v > — 1 there exists t 0 = t 0 (u) > 0 and constant > 0 such 
that 


Pi ] {t,x,y) > Cl V) ( 1 A 




(l-x)(l-y) 


(.xy + ty+F 2 y/i 


exp 


(x - yf 

2 1 


(3.1) 


whenever x, y G (0,1) and t <t 0 . 


Proof. We generally assume that t < 1 and we begin with the case when the space 
arguments are bounded away from 0, i.e. x,y > 1/32. The constant 1/32 is chosen for 
technical reasons. Since obviously TL/ 4i1 ) < Ti we can write for any Borel set A C (0,1) 
that 

[ p { y\t,x,y)m^\dy) = E^[t < Ti; R(t) e A] > E«[f < T (x/ ^y, R(t) G A]. 

J A 

Now, applying the absolute continuity property (j2.9j) we can see that the last expression 
is equal to 


jr/- 1 / 2 ) 


t < T>/ 4 ,i)> R{t) *= ^ 4 ; 


m 

x 


v+ 1/2 


exp 


- 1/4 


ds 


RM 


(3.2) 


Notice that we can estimate the above-given integral functional on {t < TA,/ 4 n} in the 
following way 

y ds i6 y 9 , 

< -T / ds < 16 • 32 2 t < 16 • 32 2 , 


R(s) 2 x 2 
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whenever t < 1 and x > 1/32 as we have assumed. Thus we have 
[ x, y)mV\dy) > c 0 


I a Kx 


yy +1 / 2 (— 1 / 2 ) /, x , 


and consequently, by (12.131) . we obtain 

(y)u O (, A 3 x(y~x/ 4 ) 

Pi ( t,X,y) > - -v mrr I 1 A 


(; xy ) u+1 / 2 


4 1 


1 A 


(1 -x)(l -y)\ 1 
t ) Vt 


exp 


\ x - y\' 

2 1 


for every t < 1 and x,y > 1/32. Observe also that since y > x > 1/32 and t is bounded 
we have 

3x(y — x/4) 1 

1 A ——-— «1« 


4 1 


(t + xy ) u+l / 2 ’ 

which ends the proof in this case. 

Now we assume that x and y are bounded away from 1, i.e. x,y < 1/4. Recall that 
the subtrahend in the Hunt formula (12.101) is given by 

ft 


r[ u \t,x,y)= I pW(t-s, l,y)q^l(ds). 




( 0 / 


' o 

Using (12.71) we obtain that for v > —1/2 we have 

ft i 






< 


> o 


(f — s + yY +l / 2 \Jt — s 


exp 


{t-s) 


H-l 


exp 


9 


32 {t - s) 


(1 -V? 
2 [t - s) 

<lxl(ds). 


Qx,l(ds) 


The function exp(—9/(32f)) is increasing on (0, U], where t\ = 9/(32(z/ +1)). Thus, 
for t <t\ we have 

ct 


ri\t,x,y) < -ryexp I ) / g^(ds). 


^+i 


9 


32t 




Estimating the last integral simply by 1 and using (12.71) together with the fact that under 
our assumptions on x and y we have \x — y\ < 1/4, we can write 

_9_ + \x-y\ 2 ' 

32 1 


r [v) (t,x,y) 

pM(t,x,y) 


< c 2 (l + xy/t) u+l/ 2 exp ( xttt + 


2t 


< c 3 


^+1/2 


exp — 


32 1 


Notice that for v G (—1, —1/2] and t < 9/16 in a similar way we can arrive at 


ri\t,x,y) 


, Y (t-s + y )-''- L / 2 

7o - S 


exp 


(i-y ) 1 

2(f - s) 




< (13/16) -i/ - 1/2 -=exp f- — 


Vt 


32 1 


< c 4 (f + a;//)^ 172 exp ( + 


9 

32f 


f- ® ). 


2t 


32t 


and consequently 
p^(t,x,y) 

Thus, in both cases (v greater or smaller than —1/2), for t sufficiently small we have 
Pi(t,x,y) =p {v) (t,x,y) (l - r -rr > YP (u \t : x,y), x,y< 1/4. 

Since (1 A (1 — x)(l — y)/t) ~ 1 in this case, the usage of (12.71) gives the result. 
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Finally, we take x < 1/8 and y > 1/4. By the Chapmann-Kolmogorov equation we 
can write 

fl/4 

Pi\t,x,y)> / p[ u \t/8,x,z)p[ u \7t/8,z,y)7n (u \dz). 

J 1/32 

Notice that we can use the previously obtained estimates since x, z < 1/4 and z,y > 1/32 
in the integral above. Consequently, 

(x-y) 2 ' 


Pi\t,x,y) > j exp 


c 6 


2 1 


where 


K(t,x,y) = 


"1/4 


' 1/32 {t + XzY+ 1 / 2 


(1 — z)(l — y) 

1 A -- j -— J exp 


K{t,x,y), 

(lx + y — 8z) 2 


14/ 


nY v \dz). 


Note that for 1/32 < z < 1/4, x < 1/8 and y > 1/4 we have 

t + xz « t + x?/, and (1 — z)(l — y)/t « (1 — x)(l — t/)/i. 
Moreover, since z > 1/32 we have 

> \ 

wY\dz 


1//4 ^ (To? + 2/ — 8^r) 2 \ . i/ f 1/M / (7t + y — 8z) 2 


exp 


'1/32 


14t 


exp 


'1/32 




2 — 7 x—y 

8 Vt 

-4,(7x+y) 
32 Vt 


exp 


14/ 

64w 2 

U~ 


dz 


dw, 


where the last equality is obtained by substituting 8z — y — lx = 8 w\ft. Finally note that 
1/4 < 7 x + y < 15/8 and consequently the lower bound of integration is non-positive 
and for t < 1 the upper bound is greater than 1/64, which implies that the last integral 
given above can be estimated from below by a constant. Thus, combining all together 
we obtain that 


K(t,x,y) < 


c 7 Vt 


(■t + xy) v+ C 2 


1 A 


( l-x){l-y ) 


which ends the proof. 


□ 


Proposition 2. For every v > — 1 we can find > 0 and Z 0 = t 0 (is) > 0 such that 


Pi\t,x,y) < 0^ ( 1 A 




(l-x){l-y) 


(xy + t) u+1 / 2 y/t 


exp 


{x - yf 
2 / 


for every x,y G (0,1) and t <t 0 . 


Proof. The Hunt formula (12.101) together with (12.71) immediately imply that 

Pi\t,x,y) < p {u \t,x,y ) 

which gives the result in the case (1—x)(l— y) > t. Thus, from now on we will assume that 
(1 — t)(1 — y) < t < 1/4. Notice that it implies that 1 — y < >ft <1/2 and consequently 
y >1/2. We split the proof into to parts. First, we assume that x < 1/2. The proof in 
this case is based on an idea to mimic the refection principle, which is true for Brownian 
motion and obviously is not for Bessel processes. However, it leads to correct estimates. 
More precisely, for every Borel set A C (0,1) we write E x[t < 7\; R(t) E A] as 

(E<">[B(i) € A}- EM[fl(t) £2-A})+ (E/)[R(() e 2 - A] - E<f>[i > T i; R(t) e ^]) 


1 1 / (x-y) 2 / 
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and denote the first part as Iii(x,t,A) and the other as K 2 (x,t, A). Note that in the 
classical Brownian motion case, the second part vanishes and the reflection principle just 
gives the formula for the transition density function of the process killed at 1. Here we 
obviously have to deal with both parts and consequently, we begin with the following 
estimates for y > 1/2 and t < 1/4 


Jv)/ J. _ o .A /o „.\ 2u+l /<-> \ v+1 


p {v \t,x,2 -y) [2 -y 


p(A(t,x,y) V V 


= 1 


2 -y\ I v {x{2-y)/t) 


y 


Iv{xy/t) 


exp 


2 ( 1 -!/) 


t 


< 1 — exp 


2(1 ~y)\, i -y 

< c r 


t J t 

Here we have just simply used the monotonicity of z u+1 and I v (z). Thus, we have 
Ki(x,t,A) = j (p {u \t,x,y)-p {u \t,x,2-y)(^—^j \m {v \dy) 


(3.3) 


< Ci 


i -y 


p( v \t, x, y) : m yV \dy). 


Moreover, for y > 1/2 and t < 1/4 we can write 


2-y\ 2v+1 P {v \t^2-y) 

y ) p (l/) (M,?/) 


l = 


y 


2-yX +1 I v {(2-y)/t) 


h(y/t) 


exp 


2(1 - V) 

t 


Now applying (I2.4p (note that (2 — y)/t > y/t > 2) we arrive at 

2v+l 

- 1 < c 2 (l - y). 


2 — y\ 2u+1 P^ v \t, 1,2 2 /) _ i < (2-y 


y j V y 

Since A C (0,1), the condition R{t) G 2 — A implies that t > Tj and consequently 
K 2 (x, t, A) = EM[t > T i; R(t) G 2 - A\ - EM[t > Tp R(t) G A]. 

The Strong Markov property gives 

K 2 (x, t, A) = E > Tp Hi(t — Ti, A)\, 

where 

#i(s, A) = E^ [R(s) 62-4]- E^ [R(s) G A] 

= j p M (s, l,y)m {v \dy) - [ p (,/) (s, 1, y)m {v \dy) 

2u+l 


1. 


(3.4) 


' 2—A 


2 -y 

y 


P {u) (t, 1,2 -y) -p H (f, 1,2/) m {u \dy). 


The estimate provided in (13.41) enable us to write 

Hi(s, A) < c 2 [ (l-y)pM(s,l,y)mM{dy) = c 2 E ( " ) [R(s)eA-,l-R{s)}. 
J A 

Now we can apply the Strong Markov property again to come back to 

1-2/ 


K 2 (x, t, A) < c 2 / (1 -y)p (v) (t,x,y)m {v \dy) <c 2 / —-^p {u \t,x,y)rn (u \dy). 
J A J A 1 

Combining all together leads to 

P ( i\t,x,y) <c 3 ^-^-p (u \t,x,y), x < 1/2, y > 1/2, t < 1/4, 
which ends the proof in this case. 


(3.5) 
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Now we assume that x G (1/2,1) and write for every Borel set A C (0,1) 

E <r>[i < Tu R(t) e A] = E <”)[* < T W1)1) ; fl(() eA]+ eM[T w1|1) < t < Tr,R(t) € A], 

(3.6) 


Intuitively, the first part should be larger then the other one, since for x,y > 1/2 we 
should expect that there are more trajectories going from x to y which do not go below 
x/A, then those visiting level x/A before reaching y. Indeed, note that the first term 
can be estimated in the same way as in the proof of Proposition [U i.e. we can apply 
the absolute continuity property, which together with the boundedness of the integral 
functional (for x > 1/2 and t < 1) appearing in (13.21) imply 

E 44 [t < T {x/Ajl y, R(t) E A] <c 4 ' p[~ 1/2) (t, x,y) dy. (3.7) 

Now it is enough to show that the second term in the right-hand side of (13.61) is signifi¬ 
cantly smaller then the expression on the left-hand side. To see this, we apply the Strong 
Markov property to write it as 

E 44 \T{ x /a,i) < Ti, H 2 (t — T( x /a,i)i x/A, A)] , 

where 


H 2 {u,x/A,A) = E % w)) [u < Tj; R(u) E A] — E«[u < T i; R(u) G A]. 

The last equality follows from the fact that R{T^ X / A ^) = x/A on {T( x /a, i) < t < Xj}. 
Since x/A < 1/2 we can apply the result given in (13.51) to estimate H 2 (u,x/A, A) in the 
following way 

E^Ju < Tp R(u) EA}<c 3 J ^1 A ) ^ ex P M^/ 4 _ V ?/( 2m )) ™/ v \dy) 

and let us denote the right-hand side above as H 3 {u,x/A, A). Notice that the function 
u -A u~ a exp(— b 2 / (2u)) is increasing for u < b 2 /(2a) and we have \x/A — y\> 1/4, thus 
u —> H 3 {u,x/A, A) is increasing for u < 1/48. Moreover, we have 

(:v - ^/ 4 ) 2 - (y - x/2) 2 > ^( 2 y - zx/a) > ^ x, y > 1 / 2 , 

and consequently, for u < t we obtain 


H 3 (u,x/A, A) < exp 



H 3 (u,x/2,A). 


Collecting all together with the fact that T^ x /2,i) < Tt x /4,i) we can write for t < 1/48 that 
E 44 [T( x /4,1) < t < Tp H 2 (t — T( x / 4>1 ), x/A, A)] 

— ^E 44 \T( x / 2,1) < t < Ti\H 3 (t — T( x /2,i), x/A, A)] 

< c 3 exp (—1/(64f)) E 44 [T( x/ 2 ,i) < t < Tp, H 3 (t - T (x/2A) , x/2, A)\ 

< c 5 exp(-l/(64t))E 4/) [T (x/ 2 ) i) < t < Tp H 2 (t - T {x / 2 ,i), x/2, A)] 

= c 5 exp (-1/(641)) E 44 [T( x / 2 ,i) < t < Tp R(t) G A] , 

the last two lines follows from Proposition [Hand the Strong Markov property respectively. 
Since obviously 

E 44 [r (x/2)1) <t<T 1 ; R(t) G 4 ] < EM [t < Tp X(f) G A] 
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we can choose t 0 < 1/48 in such a way that for t < t 0 we have csexp (—l/(64t)) <1/2 
and consequently 

E<r>[i < T,; R(t) eA] = E«[( < T MV)] R(t) ei] + <t<T 1] R(t ) e A] 

< EM[« < T W 4,1); R(t ) 6 A] + < 71; fi(«) 6 A] 

and the estimates given in (13.711 ends the proof. □ 


Theorem |Tj is now the consequence of the following. It is known that there exists 
C^ > 1 such that 

-^y(l -x)(l -r/)exp(-A^t) < G u t (x,y) < C^\l - a:)(l - y) exp(-A^f), (3.8) 

for every i,i/6 (0,1), where the upper bounds holds for t > T 0 for arbitrary T 0 > 0, but 
in the lower bounds we have only the existence of such To (see for example Theorem 3.7 in 
my Thus, the upper bounds in (11.211 are just the consequence of (13. 8 H and Proposition 
El To finish the proof it is enough to show that the result of Proposition [Tj is true for 
t < T 0 for arbitrary T 0 > 0. Let to > 0 be as in the theses of Proposition [Tj and note that 
m reads as 

Pi\t,x,y) > ci (t 0 , i/) (1 -x)(l -y), x,y G (0,1), 


whenever t G [t 0 /2,t 0 ]. Thus, applying the Chapmann-Kolmogorov equation we have 


'Pi ] (21, x,y) = I p\ u> (t, x, z)p[ u> (t, z, y)m [u) ( dz ) 




M/ 




>( Cl (t 0 ,u)Y(l-x)(l-y) / (1-zYz 


2 „2u+l 


dz 


> C 3 (t 0 ,v) 1 A 


(l-x){l-y) 


: exp 


(x - yf 

2 1 


(.xy + t ) u+l / 2 y/t 

whenever t G [to/2, to], which enable us to replace the condition t < t 0 in (13.11) by t < 2to 
and in consequence by t < Tq for arbitrary Tq > 0. This ends the proof of Theorem [I] 
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